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Abstract 

Thermodynamic  data  have  been  used  to  predict  the  dependence  of  methane  conversion  on  temperature  and  oxygen  partial  pressure  in 
mixed  conducting  membrane  reactors,  and  the  corresponding  fractions  of  water  vapor,  Ht,  CO  and  CO2.  The  relations  between  methane 
conversion,  gas  composition  and  oxygen  partial  pressure  were  also  used  to  formulate  the  oxygen  balance  in  mixed  conducting  membrane 
reactors,  with  tubular  reactor  and  continuous  stirred  tank  reactor  (CSTR)  configurations.  A  single  dimensionless  parameter  accounts  for 
the  combined  effects  of  geometric  parameters  of  the  membrane  reactor,  the  permeability  of  the  membrane  material,  and  flow  rate  at  the 
entry  of  the  reactor.  Selected  examples  were  calculated  to  illustrate  the  effects  of  steam  to  methane  and  inert  to  methane  ratios  in  the  gas 
entering  the  reactor.  The  values  of  oxygen  partial  pressure  required  to  attain  the  highest  yield  of  CO  and  H2  were  also  used  to  estimate 
the  stability  requirements  to  be  met  by  mixed  conducting  membrane  materials.  Suitable  membrane  designs  might  be  needed  to  bridge  the 
difference  between  the  conditions  inside  the  reactors  and  the  stability  limits  of  known  mixed  conductors. 
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1.  Introduction 

Methane  is  the  main  component  of  natural  gas  and  biogas. 
It  is  thus  a  very  important  energy  source  and  it  might  also 
be  converted  to  a  variety  of  useful  chemicals.  However,  the 
efficiency  of  the  actual  energy  conversion  methods  is  rela¬ 
tively  low,  and  excessive  production  of  CO2  is  contributing 
to  global  warming.  New  energy  conversion  systems  (e.g. 
fuel  cells)  are  thus  needed.  Unfortunately,  these  alternative 
systems  are  still  ill-suited  to  allow  the  direct  conversion  of 
hydrocarbons  in  fuel  cells,  thus  requiring  previous  reform¬ 
ing,  or  conversion  to  syngas,  i.e.  a  mixture  of  H2  and  CO. 

Methane  can  be  converted  by  partial  oxidation  when 
oxygen  is  provided  by  transport  through  a  mixed  conduct¬ 
ing  membrane  [1-9],  or  through  an  ionic  conductor  with 
electronic  short-circuiting  or  assisted  by  electrochemical 
pumping  [10].  The  mixed  conducting  membrane  technol¬ 
ogy  should  integrate  controlled  oxygen  separation,  steam- 
or  CCU-reforming,  and  partial  oxidation  in  a  single  step  [1]. 
The  technology  has  also  a  potential  for  energy  and  cost  re¬ 
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duction,  as  well  as  reduction  of  NOA  emissions.  Thus,  new 
mixed  conducting  materials  are  still  being  studied  [11-19]. 

Under  suitable  operating  conditions  one  seeks  the  fol¬ 
lowing  reaction; 

CH4  +  \02  ^  CO  +  2H2  (1) 

However,  equilibrium  is  unlikely  to  be  attained  easily,  and 
one  may  need  to  resort  to  somewhat  excessive  oxygen  partial 
pressure  to  attain  high  conversion  of  methane,  thus  running 
risks  of  further  oxidation  of  the  components  of  the  syngas 
mixture  to  water  vapor  and/or  carbon  dioxide,  as  follows; 

CO  +  \02  <*•  C02  (2) 

H2  +  ±02  H20  (3) 

Under  real  conditions  one  should  also  consider  the  effects  of 
methane  cracking  CH4  —*■  C+2H2,  especially  on  contacting 
metallic  surfaces  (e.g.  Ni  or  Fe);  this  shortcoming  is  one 
of  the  main  limitations  of  the  Ni/YSZ  cermet  anodes  of 
solid  oxide  fuel  cells,  and  explains  the  search  for  alternative 
SOFC  anode  materials  for  direct  conversion  of  hydrocarbons 
[20-22].  Copper  is  one  of  the  rare  exceptions  of  metals 
with  low  activity  towards  methane  cracking,  and  has  been 
proposed  for  some  of  these  concepts  [22]. 
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The  main  complications  resulting  from  methane  cracking 
in  electrochemical  systems  is  carbon  deposition,  thus  block¬ 
ing  the  methane  conversion  systems,  and  indeed  causing  a 
decrease  in  efficiency,  because  CO  can  be  valuable  as  an  im¬ 
portant  contribution  of  the  total  energy  supply,  and/or  for  the 
production  of  chemicals.  Cracking  may  be  somewhat  mini¬ 
mized  by  avoiding  metallic  surfaces,  i.e.  by  using  all  ceramic 
systems,  such  as  mixed  conducting  ceramic  membranes,  and 
by  ensuring  conditions  when  partial  oxidation  and  steam 
reforming  of  methane  are  more  favorable.  To  minimize  the 
risks  of  methane  cracking  one  should  thus  avoid  the  use 
of  metallic  surfaces  at  high  temperatures,  when  cracking 
is  thermodynamically  feasible.  Ceramic  mixed  conducting 
membranes  are  expected  to  be  more  promising  than  ionic 
conductors  with  metallic  electrodes  for  electrochemical 
pumping  or  short-circuited.  One  should  also  attempt  to  en¬ 
hance  the  kinetics  of  reactions  of  partial  oxidation  and  steam 
reforming. 

Though  the  highest  permeability  was  reported  for  per- 
ovskite  materials,  such  as  (La,Sr)(Co,Fe)03_^  [1,2],  or 
related  structures  (e.g.  La2Ni04)  [17-19],  these  materials 
show  some  critical  limitations,  mainly  in  what  concerns 
the  limited  chemical  stability  and  decrease  in  permeability 
under  reducing  conditions,  significant  lattice  expansion  un¬ 
der  high  gradients  of  oxygen  chemical  potential,  excessive 
thermal  expansion,  and  possibly  also  reaction  with  CO2. 
LaMnC>3-  and  LaCrOC>3-based  materials  possess  higher  sta¬ 
bility  but  their  permeability  is  poor.  Mixed  titanate-ferrate 
compositions  [9,11-13]  are  among  the  most  promising  ma¬ 
terials  to  meet  the  requirements  of  permeability,  stability 
and  low  price. 

A  detailed  knowledge  of  the  conditions  required  for  opera¬ 
tion  of  mixed  conducting  membrane  reactions  is  thus  needed 
for  the  selection  of  materials,  reactor  configurations,  and  to 
optimize  the  operating  conditions  in  terms  of  methane  con¬ 
version,  selectivity  of  reaction  products,  and  stability  of  the 
mixed  conducting  membranes. 

2.  Gas  composition  assuming  equilibrium  without 
cracking 

One  may  thus  estimate  the  changes  in  gas  composition  for 
this  relatively  complex  system,  as  a  function  of  temperature 
and  oxygen  partial  pressure;  this  will  indicate  the  expected 
working  conditions  required  to  reach  high  methane  conver¬ 
sion  without  the  formation  of  unwanted  products  (CO2  or 
H2O).  An  ideal  condition  may  be  calculated  for  equilibrium 
conditions  in  the  gas  phase,  as  follows: 

pCO(pH2)2  =  Kl{(p02)l/2pCn4}  (4) 

pC02  =  K2{(p02)l/2pC0}  (5) 

pH20  =  ^3{(p02)1/2pH2]  (6) 

where  pi  denotes  the  partial  pressure  of  species  (atm).  The 
Gibbs  free  energies  AGi,  A G2,  and  AG3  were  calculated 


from  standard  thermodynamic  data,  and  used  to  obtain 
the  corresponding  values  for  equilibrium  constants  K\ ,  K2 
and  K3. 

The  gas  composition  will  thus  vary  with  the  oxygen  partial 
pressure  (at  constant  temperature),  and  one  still  needs  two 
additional  conditions  to  estimate  the  gas  composition  for  a 
given  value  of  oxygen  partial  pressure.  Such  conditions  can 
be  obtained  by  taking  into  account  the  H:C  =  4  ratio  in 
methane,  or 

PH2  +  PH20  =  2(pCO  +  pCOo).  (7) 

In  addition,  the  sum  of  the  fractions  of  the  different  gases 
should  be  unit,  i.e. 

PH2  +  pH20  +  pCO  +  pCXh  +  pCU4  =  1  (8) 

A  solution  for  this  problem  can  be  obtained  on  combining 
Eqs.  (5)-(8)  to  express  the  fractions  of  CO,  H2,  CO2,  and 
H2O  as  functions  of  the  fraction  of  methane  and  partial 


pressure  of  oxygen  as  follows: 

(1/3)(1  -  PC H4) 

P  1  +  K2(p02)'f2 

(9) 

(2/3)  (1  -  /2CH4) 

P  2  1  +  K2(p02)'/2 

GO) 

(1/3)(1  —  PCH4) 

P  2  1  +  [K2(p02)l/2]~' 

(11) 

(2/3) ( 1  -  PCH4) 

p  2  \  +  [K^po2m-x 

(12) 

Substitution  of  Eqs.  (9)  and  (10)  in  Eq.  (5)  thus  yields  the 
dependence  of  partial  pressure  of  methane  on  the  oxygen 
partial  pressure: 


U  ^4)  =  f^i^02)1/2[l  +  ^2(p02)1/2] 

pCt±4 

x[l  +  /f3(p02)1/2]2  (13) 

The  dependence  of  p02  on  [CH4]  was  solved  numerically 
(for  given  values  of  pCH4),  and  the  corresponding  solution 
was  replaced  in  Eqs.  (9)— (12)  to  obtain  the  partial  pressures 
of  the  remaining  components  of  the  gas  mixture.  Fig.  1 
shows  calculations  for  a  typical  temperature  T  —  1173K. 


3.  Combined  partial  oxidation  and  water  vapor 
reforming  of  methane 

One  may  also  attempt  to  design  membranes  for  a  com¬ 
bination  of  partial  oxidation  and  water  vapor  reforming  in 
mixed  conducting  membrane  reactors.  Such  an  approach 
should  have  advantages  in  what  concerns  the  decrease  in 
oxygen  flux  requirements,  increasing  the  H2:CO  ratio,  and 
also  in  terms  of  rendering  the  gas  input  less  reducing,  and 
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Fig.  1 .  Equilibrium  gas  composition  and  methane  conversion  as  a  function 
of  oxygen  partial  pressure,  at  1100K,  for  steam:methane  ratio  0.02. 


thus  less  aggressive  towards  the  mixed  conducting  mem¬ 
brane  materials.  Under  equilibrium,  steam  reforming  should 
be  given  by 

CH4  +  H20  CO  +  3H2  (14) 

which  corresponds  to  a  combination  of  reactions  1  and  3. 

Thus,  one  might  still  combine  Eqs.  (4)-(6)  and  (8).  How¬ 
ever,  C:H  ratio  will  be  dependent  on  the  water  vapor  content, 
which  invalidates  Eq.  (7).  Alternatively,  the  quantities  of  the 
different  gas  species  can  be  expressed  as  a  function  of  the 
fraction  converted.  One  may  thus  derive  those  relations  for 
1  mole  of  methane  in  the  input  with  the  corresponding  val¬ 
ues  of  steam  to  methane  ratio  wq  =  H20:CH4,  and  an  inert 
to  methane  ratio  z.  On  converting  a  fraction  a  of  methane 
one  obtains  1  —  a  moles  of  methane,  a  moles  of  CO  +  C02, 
wo  +  2a  moles  of  H2  +  H20,  and  the  original  quantity  z  of 
inert.  From  these  quantities  one  easily  obtains  the  following 
relations  for  the  gas  composition: 


pCH4  = 

1  +  z  +  wo  +  2a 

(15) 

pCO + Pco2 = 

a 

06) 

1  +  z  +  wo  +  2a 

pR2  +  pH20  = 

wo  +  2  a 

(17) 

1  +  z  +  wq  +  2a 

On  combining  Eq.  (5)  with  Eq.  (16),  and  Eq.  (6)  with 
Eq.  (17)  one  thus  obtains: 


m  “/(l+^  +  w0  +  2  a) 

p  1  +  A2(p02)'/2 

(18) 

(wo  +  2a)/(l  +  z  +  wo  +  2a) 
p  2~  i  +  K3(Po2)'n- 

(19) 

a/(\  +  z  +  w0  +  2a) 

P  2  ~  1  +  [K2(p02)i/2]— 1 

(20) 

„  „  (u>0  +  2a)/(l  +z  +  wo  +  2a) 

,Ml°  =  i  +  [Xj(^2,./2]-.  <21) 

On  replacing  Eqs.  (14),  (18)  and  (19)  in  Eq.  (4)  one  finds 
the  additional  solution  required  to  obtain  the  value  of  p02 
for  the  assumed  fraction  reacted 

(wo  +  2  a)2a 

(1  +  z  +  wo  +  2a)2(l  -  a) 

=  Ki(p02)l/2[1  +  K2(p02)l/2][1  +  K3(P02)^2]2 

(22) 


A  Fibonacci  method  [23]  was  then  used  to  obtain  the  solu¬ 
tions  of  fraction  reacted  a  for  pre-selected  values  of  oxygen 
partial  pressure.  The  values  of  a  and  p02  were  then  in¬ 
serted  in  Eqs.  (18)— (21)  to  obtain  the  gas  composition.  The 
examples  presented  below  were  computed  using  the  follow¬ 
ing  equilibrium  constants  Kj  =  exp[— AG,/(RT)]  extracted 
from  thermodynamic  data  (e.g.  [24]) 


A4  =  14.6exp 


5.12  x  103  +  6.171og(7) 


K2  =  2.92  x  10  5  exp 


K2  =  1.37  x  10  3  exp 


T 

3.40  x  104 
T 

2.96  x  104 


(23) 

(24) 

(25) 


Typical  examples  are  shown  in  Figs.  1  and  2,  for  methane 
with  relatively  low  humidity  (wo  =  0.02),  and  for  a  mix¬ 
ture  with  water  to  methane  ratio  wo  =  0.25.  The  predictions 
show  a  slight  enrichment  of  hydrogen  in  the  syngas,  which 
may  contribute  to  increase  the  efficiency.  In  addition.  Fig.  3 
shows  a  shift  towards  slightly  less  reducing  conditions,  for 
identical  values  of  methane  conversion,  which  corresponds 
to  somewhat  less  severe  working  conditions.  However,  this 


Fig.  2.  Equilibrium  gas  composition  and  methane  conversion  as  a  function 
of  oxygen  partial  pressure,  at  1100K,  for  stearmmethane  ratio  0.25.  The 
dashed  line  represents  the  fraction  converted,  and  the  symbols  represent 
the  mole  fractions  of  H2  (H),  CO  (A),  CO2  (A),  and  H2O  (□). 
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log(Po2/atm) 

Fig.  3.  Dependence  of  the  fractions  of  FF  (B,  □)  and  CO  (A.  A)  vs. 
oxygen  partial  at  1200  K,  with  z  =  0,  and  steam:methane  ratio  wo  =  0.01 
(closed  symbols)  and  0.5  (open  symbols). 


effect  is  negligible  for  nearly  complete  methane  conver¬ 
sion.  Though  the  fractions  of  different  gas  species  are  not 
shown  in  Fig.  3,  these  can  be  calculated  by  taking  into  ac¬ 
count  Eq.  (1 8)— (2 1),  with  the  relevant  equilibrium  constants 
(Eq.  (23)-(25)).  Similar  calculations  are  shown  in  Fig.  4,  for 
methane  diluted  with  an  inert  gas,  and  wet. 


4.  Ideal  solutions  for  tubular  reactors 

In  the  concept  of  mixed  conducting  membrane  reactors 
the  gas  composition  will  be  dependent  on  the  permeability 
of  those  membranes,  and  reactor  configuration.  For  contin- 


log(Po2/atm) 

Fig.  4.  Dependence  of  the  fraction  of  FF  (B.  □)  and  CO  (A.  A)  as  a 
function  of  oxygen  partial  pressure,  at  1 1 00  K,  for  a  stearmmethane  ratio 
(too  =  0.2),  and  for  inert: methane  ratio  2  =  0  (A,  □)  and  2  (A.  B). 


Fig.  5.  Schematic  representation  of  a  tubular  mixed  conducting  membrane 
reaction. 


uous  reactors  one  might  assume  typical  chemical  engineer¬ 
ing  concepts  such  as  the  continuous  tubular  reactor,  and  the 
continuous  stirred  reactor.  In  this  section  one  considers  a 
tubular  concept,  with  a  mixed  conducting  membrane  of  in¬ 
ternal  diameter  D,  and  thickness  h  <£  D  (Fig.  5).  The  oxy¬ 
gen  permeation  rate  (moles  per  unit  area  and  per  unit  time) 
across  the  walls  of  mixed  conducting  membrane  is 

JO2  ln{(p02)f/ p02]  ^ 

J  O2  = - - -  (26) 

h 


where  JO2  is  the  specific  oxygen  permeability  of  the  selected 
membrane  material  and  h  is  the  membrane  wall  thickness. 
The  oxygen  partial  pressure  in  the  feed  side  of  the  membrane 
where  (p02)f  is  typically  0.21  atm,  in  air. 

A  molar  balance  to  oxygen  may  be  expressed  on  selecting 
an  element  of  length  SL ,  and  writing  the  input  and  output  of 
all  oxidized  species  as  follows: 


2 


JO2  ln[(pQ2)f/j?Q2] 
h 


jtDSL  =  CtQySYq 


(27) 


where  Cj  is  total  gas  concentration,  Qy  the  input  volume 
flow  rate  entering  the  tubular  reactor,  and 

Y  =  (pH20  +  pCO  +  2pC02)  (28) 

This  is  a  weighted  contribution  of  all  oxygen  containing 
species  after  accounting  for  the  expected  volume  expansion, 
i.e.  Q v  >  Qy.  Though  oxidized  species  might  also  diffuse 
up-stream,  and  reduced  species  down-stream,  one  assumed 
that  these  effects  are  negligible  when  compared  to  the  oxi¬ 
dation  process,  and  to  the  flow  rate  in  the  longitudinal  direc¬ 
tion.  This  may  not  be  the  case  when  diffusion  is  relatively 
fast,  and/or  the  oxidation  is  slow,  thus  requiring  low  flow 
rates. 

Volume  expansion  can  be  described  by  recalling  the  rele¬ 
vant  relations  between  the  contents  of  different  gas  species 
and  the  fraction  converted,  shown  above.  These  relations 
yield 

Qy  _  1  +  2a  +  z  +  wo 

Q%~  l+z  +  wo  {  ) 

and  on  replacing  this  in  Eqs.  (18),  (20),  (21)  and  (29)  in 
Eq.  (28)  one  obtains 


Y  =  (1  +  z  +  wo)_1{(wo  +  2a)  (1  +  K3~l  P02~l/2r1 

+  a  +  a(l  +  K2-lp02~1'2)-1} 


(30) 
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Fig.  6.  Simulation  of  the  dependence  of  methane  conversion  ( — ),  gas 
composition  (symbols),  and  oxygen  partial  pressure  ( — )  as  a  function 
of  dimensionless  permeability,  for  a  tubular  reactor  concept.  The  symbols 
represent  the  mole  fractions  of  H2  (H),  CO  (A),  CO2  (A),  and  H2O  (□). 


Therefore,  the  values  of  Y  can  be  related  to  previously  com¬ 
puted  sets  of  results  for  the  fraction  of  methane  converted, 
and  the  corresponding  values  of  oxygen  partial  pressure. 
Thus,  Eq.  (27)  can  be  rearranged  as  follows: 


tcDL  JO2  ^  f  Y,m' .  (p02)i  ~ 

hQyCj  JviB  L  PO2  J 


dr 


(31) 


and  solved  to  obtain  the  conditions  required  to  attain  the 
desired  conversion.  The  value  of  Y  at  the  entry  of  the  tubular 
reactor  reduces  to 
wo 


Tin  = 


1  +  z  +  W0 


(32) 


Actually,  these  calculations  were  based  on  previous  calcu¬ 
lations  of  the  dependence  of  fraction  reacted  and  gas  on 
oxygen  partial  pressure,  for  a  series  of  values  of  pC>2,  thus 
yielding  sets  of  (pC>2)k,  oik,  k  =  0,  1,  2,  3, ... ,  and  the  cor¬ 
responding  values  of  Y A  Runge-Kutta  method  [25]  was 
then  used  to  perform  the  integration  of  Eq.  (30). 

An  example  of  these  calculations  is  shown  in  Fig.  6,  for 
wet  methane  with  wo  =  0.2.  These  results  show  the  depen¬ 
dence  of  gas  composition  on  the  dimensionless  permeability 

r  =  (33) 

m:i  e°  1‘Crai 


A  being  the  exchange  area  of  the  tubular  reactor.  The  di¬ 
mensionless  permeability  accounts  for  the  combined  effects 
of  the  membrane  permeability,  thickness,  other  geometrical 
characteristics  of  the  tubular  reactor  (diameter  and  length), 
and  flow  rate.  For  example,  a  membrane  material  with  in¬ 
sufficient  permeability  requires  either  thinner  membranes, 
or  lower  flow  rates  to  attain  the  required  conversion.  Other¬ 
wise,  one  needs  tubular  reactors  with  longer  length  or  larger 
diameter.  Note  also  that  the  representation  shown  in  Fig.  6, 
also  corresponds  to  the  changes  in  composition  and  oxygen 


E 

+■» 

-5 

CN 

o 

CL 

Ui 

O 


J* 


Fig.  7.  Simulation  of  the  dependence  of  methane  conversion  ( — ),  and 
oxygen  partial  pressure  ( — )  as  a  function  of  dimensionless  permeability, 
at  1100  K,  and  for  different  values  of  the  steam:methane  ratio  (wo  =  0.01, 
0.2,  0.5). 


partial  pressure  along  the  length  of  the  tubular  reactor,  for 
given  values  of  the  membrane  diameters  and  thickness,  and 
for  flow  rate. 

The  effects  of  water  vapor  partial  pressure  at  the  entry  of 
the  tubular  reactor  are  shown  in  Fig.  7,  and  effects  of  dilu¬ 
tion  of  methane  with  an  inert  gas  are  shown  in  Fig.  8.  These 
calculations  show  that  one  might  rely  on  membrane  materi¬ 
als  with  lower  permeability,  or  use  a  lower  membrane  area 
to  attain  a  given  fraction  converted  by  promoting  simulta¬ 
neous  steam  reforming. 

Eq.  (31)  may  also  be  rewritten  to  describe  the  changes 
of  oxygen  chemical  potential  at  a  generic  distance  from 
the  entry,  along  the  length  of  a  tubular  reactor,  and  for  a 
given  conditions  of  membrane  permeability,  membrane  di¬ 
ameter  and  total  length,  and  volume  flow  rate  at  the  input, 
i.e.  for  a  given  value  of  dimensionless  permeability  J*  = 


J* 


Fig.  8.  Simulation  of  the  dependence  of  methane  conversion  ( — ),  and 
oxygen  partial  pressure  ( — )  as  a  function  of  dimensionless  permeability, 
at  1100K,  for  wo  =  0.2,  and  with  different  values  of  inert:methane  ratio 
z  =  0,  1,  and  2. 
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Fig.  9.  Simulation  of  the  evolution  of  gas  composition  (symbols),  and 
oxygen  partial  pressure  ( — ),  along  the  length  of  a  tubular  membrane 
reactor  with  a  dimensionless  permeability  J*  =  0.014.  The  symbols 
represent  the  mole  fractions  of  H2  (M),  CO  (A),  CO2  (A),  and  H2O  (□). 


7TDLJ02/hQyCT.  In  this  case 


x 

L 


(P02)f 

p02 


AY' 


(34) 


A  typical  example  is  shown  in  Fig.  9,  to  emphasize  that  the 
membrane  materials  are  more  likely  to  undergo  degradation 
under  severe  reducing  conditions  near  the  input.  One  should 
thus  rely  mainly  on  steam  reforming  near  the  entry,  by  re¬ 
sorting  to  suitable  catalysts  supported  on  materials  with  en¬ 
hanced  stability. 

The  predictions  shown  in  Figs.  6-9  are  probably  unreal¬ 
istic,  at  least  at  the  entry.  In  fact,  it  is  assumed  that  steam 
reforming  occurs  instantaneously  on  entering  the  tubular  re¬ 
actor.  The  effects  of  kinetic  limitations  are  being  simulated, 
and  will  be  dealt  with  in  a  separate  report. 


Fig.  10.  Predictions  of  methane  conversion,  gas  composition  and  oxygen 
partial  pressure  on  dimensionless  permeability,  for  a  CSTR  reactor  con¬ 
cept.  The  symbols  represent  the  mole  fractions  of  H2  (H),  CO  (A),  CO2 
(A),  and  H20  (□). 


Eq.  (30)  can  be  used  to  obtain  the  value  of  r0ut  for  the 
required  conversion,  and  7;„  =  wo/(l+z  +  wo). 

Eq.  (36)  can  thus  be  used  to  describe  the  dependence  of 
fraction  crystallized  on  the  permeability  of  the  membrane 
material,  geometric  characteristics  of  the  reactor,  and  flow 
rate,  all  included  in  the  dimensionless  permeability: 


JOiA 

hCjQ0 


Thus,  the  dependence  of  fraction  reacted  on  the  permeabil¬ 
ity  is  easily  obtained  by  inserting  the  previously  computed 
dependence  of  fraction  reacted  on  the  oxygen  partial  pres¬ 
sure.  One  typical  example  is  shown  in  Fig.  10.  Though  this 
CSTR  concept  is  less  sensitive  to  a  very  reducing  input,  it 
is  very  likely  to  require  much  larger  volume  to  area  ratios, 
and  thus  lower  efficiency  for  partial  oxidation. 


5.  Model  for  CSTR  reactors 


6.  Stability  requirements  of  membrane  materials 


A  simpler  balance  can  be  derived  on  assuming  a  perfectly 
stirred  continuous  reactor,  with  volume  V,  with  a  volume 
flow  rate  Q y  at  the  entry,  and  an  area  A  of  mixed  conducting 
wall  for  the  supply  of  oxygen.  In  this  case,  the  relevant 
balance  of  oxygen  reduces  to 


2A 


JP2  In[(p02)f/p02] 
h 


CiQvipY bO  +  PCO 


+  2pC02]  -  CtGv{(/T120)o  +  (pCO)o  +  2(pCO2)0) 


(35) 


This  can  be  combined  with  Eqs.  (18),  (20)  and  (21)  and 
re-arranged  to  obtain  the  following  generic 


JQ2A 

hQyCj 


0.5  In 


(pP2)f 

p02 


(Tout 


Tin) 


(36) 


Figs.  1-4  show  that  the  highest  values  of  the  fractions 
of  H2  and  CO  are  attained  for  p02  close  to  10-21  atm, 
at  1100K.  The  corresponding  conversion  of  methane  is 
already  in  the  order  of  95%.  On  shifting  the  operating  con¬ 
ditions  to  more  oxidizing  one  may  thus  expect  the  onset 
of  completely  oxidized  species  (i.e.  C02  and  H20).  The 
conditions  required  to  attain  the  highest  yield  of  H2  and 
CO  are  slightly  temperature  dependent  (dashed-dotted  line 
in  Fig.  11),  whereas  effects  of  humidity  (Fig.  3)  and  of 
dilution  with  an  inert  gas  (Fig.  4)  are  relatively  weak. 

The  results  shown  in  Fig.  11  are  also  a  guideline  for 
assessing  the  stability  requirements  of  mixed  conducting 
membrane  materials.  Note  that  known  mixed  conductors 
have  limited  stability  under  reducing  conditions.  Typical 
examples  are  materials  with  perovskite  or  related  structure 
types,  such  as  LaCoOs  or  LaFeOs,  with  partial  substitution 
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T  (K) 


Fig.  11.  Predictions  of  conditions  required  to  obtain  the  highest  yield 
of  H2  and  CO  by  partial  oxidation  of  methane  with  steam: methane 

ratio  wo  =  0.25  ( - ),  and  typical  stability  limits  of  prospective  mixed 

conductors,  and  stability  of  other  oxides  (NiO,  CoO  and  FeO)  added  to 
mixed  conductors  based  on  lanthanum  gallate  or  titanates. 

of  Ln3+  by  Sr2+,  and  La2Ni04.  Their  stability  ranges 
were  thus  assessed  using  results  reported  in  the  literature 
[26-29],  and  are  compared  to  the  expected  working  condi¬ 
tions  required  for  the  conversion  of  methane.  These  results 
suggest  that  one  should  seek  candidate  materials  based 
on  LaFe03  and  possibly  also  LaMn03  with  suitable  com¬ 
position  changes,  to  enhance  the  transport  properties.  For 
example,  (La,Sr)Fe03_3  might  be  candidate  materials  with 
a  compromise  between  permeability  and  stability.  Other 
promising  mixed  conductors  are  based  on  lanthanum  gal¬ 
late  [14,15]  or  titanates  [11-13]  with  additions  of  Ni,  Co 
or  Fe.  In  these  cases,  the  stability  limitations  are  proba¬ 
bly  determined  mainly  by  the  reducibility  of  the  additives. 
The  stability  limits  of  single  oxides  (e.g.  NiO,  CoO  and 
FeO)  were  thus  computed  from  thermodynamic  data,  and 
are  shown  as  guidelines  for  materials  containing  these 
oxides. 


7.  Conclusions 

A  combination  of  analytical  and  numerical  methods 
can  be  used  to  describe  the  processes  of  partial  oxidation 
and  steam  reforming  of  methane,  under  equilibrium  con¬ 
ditions,  and  as  a  function  of  oxygen  partial  pressure.  This 
dependence  of  fraction  of  methane  converted  and  the  cor¬ 
responding  fractions  of  CO,  FFO  and  CO2  can  be  used 
to  derive  solutions  for  the  steady  state  regimes  in  typical 
tubular  reactor  and  CSTR  reactor  concepts  with  mixed 
conducting  walls.  A  single  dimensionless  parameter  may 
account  for  the  combined  effects  of  oxygen  permeability  of 
the  reactor  wall,  geometric  characteristics  and  volume  flow 
rate.  Selected  examples  showed  that  partial  steam  reform¬ 


ing  may  be  crucial  to  avoid  excessively  reducing  conditions 
at  the  entry,  and  to  protect  the  mixed  conducting  mem¬ 
brane  materials  exposed  to  those  conditions.  Slow  surface 
exchange  also  causes  a  gradient  of  chemical  potential  at 
the  permeate  side  of  the  membrane,  and  might  contribute 
to  bridge  the  gap  between  the  stability  limits  of  membrane 
materials  and  excessively  reducing  conditions  inside  the 
reactor. 
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